








Turbulence in Electrostatic 
Ion Acoustic Shocks 


PPG-149 


R.W. Means, F.V. Coroniti 
A.Y. Wong and R.B. White 


May, 1973 



TURBULENCE IN ELECTROSTATIC ION ACOUSTIC SHOCKS 


R. W. Means, F. V. Coroniti, A. Y. Wong, and R. B. White 
Department of Physics, University of California, Los Angeles, California 90024 • 


Abstract 

Three types of collisionless electrostatic ion acoustic shocks 
are investigated using the University of California, Los Angeles, 
double plasma (DP) device: a) laminar shocks; b) small amplitude 

turbulent shocks in which the turbulence is confined to be upstream 
of the shock potential jump; c) large amplitude turbulent shocks 
in which the wave turbulence occurs throughout the shock transition. 
The wave turbulence is generated by ions which are reflected from 
the shock potential; linear theory spatial growth increments agree 
with experimental values. The experimental relationship between 
the shock Mach number and the shock potential is shown to be incon- 
sistent with theoretical shock models which assume that the electrons 
are isothermal. Theoretical calculations which assume a trapped 
electron equation of a state and a turbulently flattened velocity 
distribution function for the reflected ions yields a Mach number 
vs potential relationship in agreement with experiment. 



I. 


INTRODUCTION 


Collisionless laminar electrostatic ion acoustic shocks were first 
predicted by Moiseev and Sagdeev.^ They showed that the Debye length, disper- 
sion limited solitary wave calculated using isothermal electron-cold ion 
fluid theory could be converted into a shock transition by the inclusion of a 
small fractional density of ions which are reflected by the shock potential. 

They suggested that these reflected ions might stimulate a counter-streaming 
instability thus introducing turbulent dissipation into the laminar shock 
structure. This theory predicted steady shock solutions only for Mach numbers 
below 1.6. (In the shock frame, the Mach number is defined as the upstream 
flow speed divided by the long wavelength linear ion acoustic speed.) For a 

Maxwellian distribution of upstream ions and isothermal electrons, Bardotti 
2 

and Segre showed that the maximum or critical Mach number decreased with 

decreasing electron to ion temperature ratio (T^/T^) , and that for T^/T^ < 12.5, 

no steady shock solutions were possible. In the cold ion limit, however, 

3 4 

Montgomery and Montgomery and Joyce demonstrated that if the shock were 
treated as an ion acoustic Bernstein-Greene-Kruskal mode,"* then the freedom 
to adjust the distribution of electrons trapped in the shock potential per- 
mitted shock-like solutions to be constructed for any Mach number. In numerical 
simulation experiments Forslund and Shonk^ and Forslund and Freidberg^ investi- 
gated high Mach number, high T g /T^ shocks and found steady laminar shocks with 
Mach numbers between 3 and 6 in which the trapped electron distribution function 
was essentially flat. The high T^/T^ ratios and one-dimensional property of the 
simulations precluded a study of reflected ion driven turbulence. As an alterna- 

g 

tive to the laminar shock theory, Tidman developed a fully turbulent shock model 
in which the interpenetration of cold upstream and shock heated downstream ions 
generated ion acoustic turbulence in the shock front. Tidman reasoned that 



-3- 


for high Mach numbers the laminar shock structure would break and turbulent 

dissipation would be required for a shock transition. However, the existence 

of high Mach number laminar shocks in the numerical simulations leaves open 

the possibility that ions reflected from a steady shock potential might be 

the dominant source of turbulence in electrostatic shocks. 

Laboratory investigations of electrostatic shocks began with Taylor et 
9 

al. who studied low Mach number shocks in high T^/T^ plasmas. They found 

steady laminar shock profiles that were qualitatively similar to the predictions 

of Moiseev and Sagdeev.^ Means and Wong^ and Stern and Decker^ extended 

these studies to higher Mach numbers and low T /T. ratios and found that a 

sufficient density of ions was reflected from the shock potential to generate 

short wavelength, high frequency ion acoustic turbulence in the shock front. 

Reflected ions have also been detected in experiments on high Mach number 

12 

magnetosonic shocks. For Alfven Mach numbers exceeding 3.0, resistivity 

alone is inadequate to provide the dissipation required by the magnetosonic 

13 

shock Rankine-Hugoniot relations. Several authors have proposed that an 

electrostatic ion acoustic subshock develops on scale lengths short compared 

with the magnetic transition length and produces the requisite additional ion 

dissipation. In magnetosonic shock experiments, however, the Debye length 

ion acoustic subshock is not directly detectable. Hence, the properties of 

electrostatic shocks and reflected ion-driven wave turbulence are best 

investigated in magnetic field-free plasmas where Debye scale lengths are 

directly accessible to ordinary probe techniques. 

This paper reports the results of turbulent electrostatic shock experiments 

14 

performed on the University of California, Los Angeles, double plasma machine. 
Section II discusses the experimental techniques and the laminar shock data. The 
experimental relationship between the Mach number and the maximum shock potential 
is compared with several models. The experimental values are found to lie above 
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those predicted by Moiseev and Sagdeev for isothermal electrons and cold ions 
and below those predicted by Forslund and Shonk^ for a trapped electron 
distribution with cold ions. The presence of warn ions has a significant 
affect on the trapped electron model, but has a minimum effect on the isothermal 
electron model. 

Section III discusses the turbulent shock structure. Three basic types 
of shocks were observed: 1) laminar shocks; 2) turbulent shocks in which the 
turbulence was small and did not significantly interact with the shock front; 

3) turbulent shocks in which the turbulence became large and interacted 
strongly with the shock front. The three-dimensional nature of the turbulence 
induced by the reflected beam is also discussed along with its growth rates 
and spectrum. 

Section IV is a theoretical discussion of the ion beam plasma stability 
problem. Two models are discussed: the Maxwellian beam and the truncated 

Maxwellian beam. Both of these models are approximations to the actual ion 
beam reflected by the shock potential. The group velocities of the unstable 
waves are also investigated. Qualitative agreement with experiment is found. 

Section V is a discussion of a model of the turbulence in the shock front. 
We show that the Mach number vs maximum shock potential relationship cannot 
be explained by a theory assuming isothermal electrons even if the maximum 
turbulent modification of the reflected ion distribution is considered. The 
experimental relationship is well approximated by assuming a maximally trapped^ 
electron distribution and a turbulently constructed plateau for the reflected 
ion distribution. Since trapped electrons are not directly detectable in the 
experiments, the theoretical agreement with experiment tends to confirm that 
trapped electrons do indeed determine the experimental shock structure. 
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II. EXPERIMENTAL TECHNIQUES AND MEASUREMENTS 

The UCLA double plasma machine consisted of two dc discharge plasmas 
produced in separated cylindrical metal walled chambers. Each chamber had a 
diameter of 30 cm and a length of 30 cm. They were also insulated from each 
other and separated by a wire grid with mesh spacing less than the Debye length. 
The grid was biased negative so that the electrons in one plasma were isolated 
from electrons in the other plasma, while the ions could flow from one chamber 
to the other. The plasmas were in good electrical contact with the walls of 
the chambers which served as the anodes of the discharge current. This fixed 
the plasma potential in each chamber near the value of the wall potential 
(usually several volts above the wall potential). In practice, the walls of 
one chamber were grounded and a signal source was connected to the walls of 
the other chamber. The electron temperature in the plasma was controlled 

14 

between 0.5 eV and 5.0 eV by means of a small second anode in each chamber. 

The temperature was measured by means of a swept Langmuir probe. The ion 

9 

temperature was measured by means of an electrostatic energy analyzer. It 
was not variable and was measured to be 0.1 eV. This value was very close to 
the resolution of the instrument. 

The plasma was a weakly ionized argon plasma with a variable density 
8 10 3 

between 10 and 10 /cm . Operating conditions were usually with the density 

9 -4 

at 10 /cm and a neutral pressure of 2.10 Torr. The ion charge exchange 

length was approximately 10 cm. 

A.. Shock Parameters 

Since the Mach number of a shock is defined as the ratio of the velocity 
of the shock to the phase velocity of the small amplitude ion wave, it was 
necessary to produce and measure small amplitude ion acoustic waves in the 
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double plasma device. The ion acoustic velocity was found to be slightly less 

1/2 

than that predicted by the linear theory c g = [ (T^ + Y^T^)/M_^] } where is 

an effective ion ratio of specific heats. This effect was more pronounced at 
low electron temperatures and is not completely understood. The discrepancy 
might be due to the finite size of the device and the non-Maxwellian distribu- 
tion function of the electrons. The device operates with a considerable 
density of primary discharge electrons in the plasma. In a larger double 
plasma device, the ion acoustic phase velocities have also been measured and 
the deviation from linear theory was not quite as large as in the smaller 

device. The Mach numbers quoted in this paper are defined as the ratios of 

the two experimentally measured velocities. 

The shocks were generated by producing a large amplitude density step at 
the interface between the two plasmas. The wall potential of the rear chamber 

was raised by means of a voltage ramp in a time comparable to several ion 

plasma periods. This forced the ions in the rear chamber to flow across the 
interface into the front chamber where a local density excess was produced. 

The excess ion charge and the electrons in the front chamber interacted self- 
consistently to change the ion flow energy into a wave front, which then 
propagated away from the interface at supersonic speeds. 

Figure 1 is an oscilloscope trace of a typical laminar electrostatic shock. 
It consists of electron density measurements as a function of time at several 
distances from the generating region. The reflected beam can clearly be seen 
in the foot structure ahead of the shock. The Mach number, density step, 
potential step, beam velocity, and beam density are all interrelated parameters 
that were experimentally measured. While the Mach number is directly measurable, 
the other quantities depend to some extent upon the theoretical model used to 


reduce the raw data. 
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Mach Number 

The experimental relationship between the Mach number and shock potential, 
ec()/T e , is plotted in Fig. 2 for a large number of observed shocks. The 
electron temperature was varied from 0.5 eV to 5.0 eV. The Mach number 
dependence on electron temperature was not discernible from the experiments. 

The scatter due to uncertainties in determination of the Mach number and 
shock potential were greater than any systematic shifts detected when the 
electron temperature was varied. 

Two theoretical curves are also shown on Fig. 2. Both are cold ion 

theories (T /T. = °°) . Moiseev and Sagdeev^ assumed isothermal electrons and 
e x 

arrived at the result 

[exp(ei|>/T ) - 1.0] 2 

- . ( 1 ) 

exp(ecf)/T ) - 1.0 - e<j>/T 
G G 

Forslund and Shonk assumed a trapped electron distribution function with the 
equation of state 


P(e(|)/T e )/P 0 = 2(e(j)/T e ) 1/2 //? + exp (ecj>/T e ) erfc [ (e<j>/T e ) 1/2 ] + 4 (e<}>/T e ) 3/2 /3/if 


and arrived at the Mach number given by 


M 2 = 1 (P - 1-0) 

2 (P - 1.0 - e<J>/T e ) 


The experimental data lie between these two theoretical curves. 

The effect of the finite electron to ion temperature ratio was investi- 

1 2 

gated by Bardotti and. Segre ’ for the case of isothermal electrons and 
Maxwellian ions. In order to compare their calculations with this experiment, 
the Mach number in their published data must be divided by the ion temperature 
correction factor for linear ion acoustic waves. This factor was computed 
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by Bardotti and Segre^ to be (1.0 + 3.4 T\/T e ) . When this was done, the 
resulting relationship between Mach number and shock potential differed only 
slightly from that predicted by Moiseev and Sagdeev for cold ions, and hence 
disagrees with experiment. 

The finite ion temperature corrections to the Mach number — potential 
relationship for the trapped electron distribution were calculated assuming 

r\j 

that the reflected ion distribution was a truncated Maxwellian. For T /T. % 

e l 

15 to 50, the theoretical curves agree roughly with the experimental data 
(Fig. 3). However, no direct observation of trapped electrons has been possible 
in the experiments. Hence, the question remains of whether turbulent relaxa- 
tion of the reflected ion distribution could sufficiently modify the isothermal 
electron shock theory to bring its predictions into agreement with experiment. 

In Sec. V.O we develop a model for the turbulent dissipation of the reflected 
ions and demonstrate that only the trapped electron distribution yields 
agreement with experiment. 


Reflected Ion Density and Velocity 

The number of ions reflected from the shock can be computed from energy 

considerations. In the wave frame, the shock front is a potential barrier of 

1/2 

height e<J>, so that all particles with velocity less than (2e<p/M^) will be 
reflected. Let us assume that the incoming ions are Maxwellian with a tempera- 
ture T^; then, the number density of reflected particles is 


*0 /rr 


dv exp(-v^/a^) 


(4) 


i v. 


1 /2 1/2 
where a. = (2T./M.) is the ion thermal speed, v 1 = [MC - (2e<j>/M.) ], 

1 IX IS 1 


V 2 = MC g , and n^ is the upstream plasma density. The average velocity of the 


reflected particles is also calculated straightforwardly to get 
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2 2 

dv vexp(-v /a^ ) 


v, E <v> = 

b v. 


dv exp(-v^/a^) 


(5) 


These relationships, (4) and (5), are indeterminate for a given shock 
wave unless the relationship between Mach number and shock potential height 
e<p/T is known. For the experimental relationship M = 1.0 +0.6 etp/T^, the 
density of reflected ions was calculated and plotted in Fig. 4 for the tempera- 
ture ratios T e /T^ = 10 and 30. The predicted beam densities from Sagdeev's 
cold ion relationship between Mach number and shock potential are also shown. 
Figure 5 shows the calculated beam velocities versus the shock potential jump 
for these temperature ratios. The beam velocity was measured by two 
independent methods. The velocity of the foot structure in front of the shock 
wave was measured directly on the oscilloscope traces of the data such as that 
in Fig. 1. It was also measured with an energy analyzer by sampling techniques. 
These measurements were generally within 15% of one another. 

The measured beam densities agree with those values predicted by the 
experimental relationships between Mach number and shock potential height . 

This demonstrates the self-consistency of the measurements, and that the beam 
ions ahead of the shock come as a result of reflection and not from other 
sources (e.g., streaming ions from the generating region). There were only 
small differences in the beam velocities predicted by the different theoretical 
models. The experimental uncertainty in beam velocity was larger than such 


differences . 
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in. TURBULENT SHOCK STRUCTURE 

Three qualitatively different shocks were observed in the course of the 
experiments: 1) laminar shocks in which no turbulence existed; 2) turbulent 

shocks in which the turbulence remained small (much less than the shock 
height) and did not significantly interact with the shock front; 3) turbulent 
shocks in which the turbulence grew to a level comparable to the shock height. 
There was considerable interaction between the turbulence and the shock front 
in the third type. 

Earlier, Fig. 1 was referred to as a typical laminar shock. It has few 
reflected ions, a small Mach number, a periodic trailing wave train, and no 
turbulence. It is the electrostatic shock discussed in most computer simula- 
tion studies and in most theories. 

Figure 6 is an oscilloscope trace of a turbulent shock in which the 
turbulence remained small and did not significantly interact with the shock 
front. The turbulence was primarily confined to the reflected beam region 
which extended ahead of the shock front. The oscilloscope trace is a super- 
position of many shocks and the turbulence is shown as the blurred region in 
the photograph. (The system was run at a repetition rate of 1 kHz.) Each 
trace is a time profile of the shock at a different distance from the interface 
between the two plasmas. At x = 2 cm, the shock appeared laminar; the reflected 
beam was evident and there existed a trailing wave train. Down the chamber at 
x = 4 cm, the first evidence of turbulence appeared in the foot structure 
ahead of the shock front. Farther down the chamber the turbulence became more 
pronounced and the trailing wave train disappeared. The turbulence always 
remained much less than the main shock front height. 

Figure 7 is a set of oscilloscope traces of the third type of shock. The 
turbulence grew to a level comparable to the shock front and interacted strongly 
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with it. Since the display is a superposition of many shocks, the turbulence 
again appears as the blurred region in the photograph. Individual shocks 
were studied on a storage oscilloscope and the turbulence became more apparent. 
These oscilloscope traces are shown in Figs. 8 and 9. Again the shock is 
laminar close to the interface region.^ The reflected beam was present and 
there was no trailing wave train. As the shock traveled down the chamber, the 
turbulence grew to a level comparable to the shock potential. It was no longer 
confined to the reflected beam region and the shock had a broad turbulent 
transition region between the upstream and downstream states. 

Figure 10 is a plot of the normalized amplitude of the shock front as a 
function of distance for the three types of observed shocks. There was a decay 
of the shock as it propagated down the chamber due to geometrical effects and 
charge exchange. The laminar shock and the small amplitude turbulence shock 
appeared similar in this respect. On the other hand, the large amplitude 
turbulence shock decayed more rapidly. As the shock propagates into the target 
plasma, it sweeps-up target plasma ions which then eventually form the 
reflected ion beam. This time-dependent energy exchange between the shock and 
reflected ions is analogous to linear Landau damping. Since the large amplitude 
turbulence shocks generally have larger shock potentials and occur for lower 
T e /T^, they reflect a larger fraction of the upstream ions than do the laminar 
and small turbulent amplitude shocks; hence, the large amplitude turbulent 
shocks should undergo a larger spatial damping. Another aspect of the turbulent 
shocks is that the turbulence appeared to be almost stationary with respect to 
the shock front. This will be discussed further in Sec. IV where the ion beam 
plasma dispersion relation is considered. 

The growth rates of the turbulence and the spectra have been measured and 
reported earlier. ^ It was possible to make detailed measurements of growth 
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rates at a particular frequency and wavelength, with the direction of propaga- 
tion known, by propagating a small ion acoustic wave packet ahead of the 
shock. Qualitative agreement was found with the ion beam-plasma theory 
developed in Sec. IV. Quantitative agreement with the theory was found only 
within a factor of two. The limitations of this theory are discussed in Sec. IV. 

The observations of turbulence in these experiments depended fundamentally 
upon the three-dimensional nature of the plasma. The reflected ion beam was 
able to interact with ambient noise ion waves which were traveling at angles 
to the beam velocity. This allowed a resonant interaction between the beam 
and the waves which resulted in the growth of the turbulence. The average 
propagation direction of the turbulent spectra in front of the shock wave was 
difficult to measure since parameters determining that direction, such as beam 
velocity, could not be controlled over a wide range of values. To surmount 
these difficulties, the double plasma machine was used to create a steady state 
ion beam with a given velocity traversing the front chamber by biasing the wall 
of the rear chamber positive. This approximated the region in front of the 
shock wave. Under such conditions there was ion acoustic turbulence in the 
front chamber. ^ The direction of wave propagation was measured by a disk 
shaped Langmuir probe with a diameter large compared to the wavelength of the 
dominant turbulent modes. A wave which propagated parallel to the face of the 
probe was averaged out over the surface. A wave which propagated perpendicular 
to the probe face was detected. Rotating the probe until the maximum signal 
was obtained gave the most probable direction of wave propagation. This angle 
is plotted in Fig. 11. Two theoretical curves are also shown in Fig. 11. 

These will be discussed further in Sec. IV. The most important point of the 
data shown is that the largest amplitude waves propagated at large angles to 
the beam. This emphasized the critical nature of the three dimensionality to 
the observation of the turbulent shocks. 
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In the one-dimensional computer simulations of electrostatic shocks by 

Sakanaka, et al.,^ Mason , ^ and Forslund and Shonk,^ the turbulence due to 

20 

reflected beam has been seen only once. We believe that this was primarily 
due to the one-dimensional character of their plasmas. Furthermore, in those 
cases where there was a possibility that turbulence could develop in a one- 
dimensional situation (i.e., at low temperature ratios), only a small number 
of computer simulations were run, and the shocks were followed for only a 
small distance. (Much of the simulation has been done at high temperature 
ratios. Only a few ions were reflected in these cases.) From the profiles of 
the turbulent shocks in Figs. 6 and 7, it is clear that the shock wave had to 
propagate a considerable distance into the plasma before the turbulence 
developed. This was simply due to the fact that the beam had to exist for a 
finite length ahead of the shock in order for the wave to grow over many wave- 
lengths. The computer simulations usually stopped long before this could 


happen. 
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IV. ION BEAM-PLASMA STABILITY ANALYSIS 


In general, the plasma ahead of the shock Is unstable to ion acoustic 

waves for an ion beam within only a limited range of parameters. For small 

beam densities and T > T., the interaction between the wave and the beam is 

e x 

primarily a resonant interaction. For instability to occur it is necessary 
for the slope of the distribution function to be. positive. In a one-dimensional 
plasma with Maxwellian distribution functions for the plasma and the beam this 
requirement is 


r(v.-c Q ) 

D S 


exp 


< -<w 


a . 

1 


2T . 
1 


> 1.0 


( 6 ) 


where is the veam velocity, C g is the ion acoustic velocity, r is the ratio 
of beam density to ion density, and a^ = [21^/M^] is the ion thermal velocity. 
Equation (6) states that must be greater than C g for an instability to occur. 
The beam velocity, however, cannot be too large, since then the term in the 
exponential becomes negative and dominates. An upper limit to the beam 
velocity is near 2.0 C g . Thus, for a one-dimensional plasma the beam velocity 
must be approximately between the ion acoustic speed and twice the ion acoustic 
speed. 

The limitation that the beam velocity must be approximately between the 
ion acoustic speed and twice the ion acoustic speed is true only in a one- 
dimensional plasma. For a three-dimensional plasma in which waves can travel 
at an angle 0 to the beam, this condition is relaxed to that of C g 'V cos0 

'V* 2C , or 
s 


cos 1 (C /V.) £ 0 £ cos (2C /V. ) 
s b s b 


(7) 


The unstable ion wave growth rates and the onset of turbulence were com- 
puted by the three-dimensional dispersion relation considered by Fried and Wong 
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for a homogeneous beam plasma system 


0) 

e(k,w) = 1 - l -y- 


t dv k • 


Sf . 
ia 


a k 


- 8v 


v - w 


= 0 


( 8 ) 


where species a has distribution function and plasma frequency If 

Maxwellian distribution functions are assumed and we specialize to the case 
where the ions of the plasma and the beam have equal mass and temperature, 
Eq. (8) becomes 



r ' 


q)-k»V, 

a) 

+ r Z' 

b 

ka . 

ka . 

l 

7 


i 


- 2(l+r) ~ 
e 


(9) 


i 

2 1/2 

where is the beam velocity, k^ is the ion Debye wave number (4im^e /l\) , 

i 

and Z'(s) is the derivative of the plasma dispersion function. It has also 
been assumed that the phase velocity u)/k is small compared with the electron 
thermal velocity. The Maxwellian assumption for the beam ions is only an 
approximation to the actual case. 

The ion waves tend to propagate in the direction which will maximize 
their growth rate. This allows us to solve the dispersion relation for k 
parallel to the beam direction as if in a one-dimensional theory. The maximum 
growth rates obtained can be projected back onto the angle at which the wave 
propagates. For example, if for a given beam density and temperature ratio 
it was found that one-dimensional growth rate was maximum near the beam velocity 
= 1.1 C g , and if the actual beam in the plasma was traveling at 2.0 C g , then 
the direction of propagation of the maximally unstable wave would be given by 
0 = cos ^(1. 1/2.0) = 57°. Figure 12 is a plot of maximum growth rate (maximized 
as a function of frequency) vs beam density at V^/C = 1.0. Since the growth 
rate was generally a maximum near V^/C^ = 1.0 these data can be also considered 
as a density stability limit graph. 
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The group velocity of the turbulent waves was experimentally found to be 
near the shock wave velocity, especially for the type 2 turbulent shocks. In 
the presence of unstable modes the group velocity is not a theoretically well 
defined quantity. In the limit of small growth rates, however, an expression 
can be arrived at from a variational principle which yields 






( 10 ) 


From Eq. (8) this gives the group velocity in the beam direction 


2 3 

2k a. 2 

— 7 ; — — cos0 + rV, sin 0 Z" 
2 b 

O) 


U)-k »V. 
— r — b 

k a . 
r i 


6oo CjO q , 

-rr = t — cos0 + 
6k k 

r 


Z" 


f 

li) 

+ rZ" 

w-kr • 4 j 

k a. 

k • a , 

r i 


r i 

k ■/ 


( 11 ) 


Two group velocity versus frequency curves are plotted in Fig. 13 for typical 
experimental situations. Linear theory with no ion beam present predicts 
group velocities on the order of 0.5 C g . Figure 13 shows that the presence of 
a small beam increases the group velocity. Hence, the observed tendency of the 
ion acoustic turbulence to propagate with group speeds comparable to the shock 
speed is qualitatively consistent with ion beam instability theory. However, 
it has not been possible to differentiate between type 2 and 3 turbulent shocks 
on the basis of the present linear stability analysis. First, the experimental 
reflected ion distribution is not Maxwellian, whereas the group speed calcula- 
tions are very sensitive to the shape of the reflected ion distribution function 


through the Z" 


co-k • V 
— r — b 

ka , 


term. Second, the reflected beam-plasma interaction 


is not completely steady with the beam, shock, and turbulence all interacting 


as the shock propagates down the chamber. Finally, additional detailed measure- 
ments of the angular distribution of the turbulent waves and the spatial 
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structure of the reflected Ion distribution function are needed to more 
clearly differentiate between the two types of turbulent shocks. 


A. Stability Analysis of a Truncated Maxwellian 

In the shock frame, if the upstream ions have a Maxwellian distribution, 
ions with kinetic energy per unit charge less than the maximum shock potential 
will be reflected back upstream. If these ions are specularly reflected from 
a steady shock front, the upstream ion distribution in the direction normal 
to the shock front (x direction) would be a truncated Maxwellian between the 
edge velocity -Vq(x) = - [2e/M^ (<J> (0) - (j>(x)) ] 5. v x H 00 given by 


f. = 


3/2 
it a . 


exp(-v 1 2 /a ± 2 ) 


/ 

< 





0 


[ 2 2e<J)(x)l 

1/2 

2 

V x + M. 

1 

" u o 

* 


-v 0 (x) < v x < 


-V < 

x — 


-v 0 (x) 


( 12 ) 


00 


<{)(x) is the local shock potential, <f>(0) is the maximum shock potential assumed 
to occur at x = 0, u^ is the upstream flow speed, directed toward positive x, 
v ± is the velocity magnitude perpendicular to the shock normal, a^ is the ion 
thermal speed and f^ is normalized to unity. In order to investigate the 
stability of f^ to ion acoustic perturbations, we assume that the. waves and 
linearized electron and ion distribution functions can be expanded in WKB 
eigenfunctions (see Sec. V) and thereby obtain the local electrostatic 
dielectric function 


exp(e<j)/T ) ^p. 

e(u),k) = 1 + — - — = — — + — —■ 


3f . , 

i ,3 


- 9v 


d v 


to - k ♦ v 
2 


(13) 


2 2 47Tn e e 

where = c /w , c = T /M . , u = — — with n being the upstream electron 

I) s s c I Pj e 
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density. Introducing the above ion distribution function, performing the 
2 

integrals over d v^, and integrating by parts with respect to dv^, the 
dispersion relation can be written as 


0 = k + exp 


e<f>(x)' 

T , 2 
e k 

00 

f j t 7 

0) - k M t 
x * 

T 

e 

v y 

T. k,k 

1 - X 

dt at 

- t- 

k, (2T./T ) 1/2 



i 

" 16 


2T . 
1 


k- G <-v z 

x 


W + k M t. 
x * 1 


kj. (2T ± /T e ) 


1/2 


(14) 


where the wave numbers k , k,, and k have been normalized to A n , the frequency 

to has been normalized to oo , and M. = u rt /c is the sonic Mach number. The 

p . * Os 

function G(t) is given by 


G(t) = 



/ 

l 

CN 

O 
0 ! 

- 

2 2e<j>(x) 

1/2 n 

- i 

2 

- < 

— exp 
& 

2 

a. 

i 

- 

t + 

M T 
I e ; 

4 


-t 1 < t < “> 


t < -t. 


(15) 


where t^ = 1/M^ [2e(<|>(0) - ^(x)}]^^ 2 . 

The dispersion relation was solved as a function of shock potential 4>(x) 
for real w (in the shock frame) and k, , complex k^, T g /T^ = 20, and = 1.46. 
Figure 14 shows the spatial growth rate Im k vs <j>(x) for Re k =1, and 

X X 

variable k, . The growth rates are largest near the maximum shock potential, 
as expected, since the relative velocity between the incident and reflected 
ions is smallest there. Furthermore, the oblique waves k^ > k^ are most 
unstable, as in the homogeneous calculations of Sec. IV-A. The very large 
growth rates indicate that the reflected ions will rapidly lose their momentum 
to the unstable ion waves and their distribution function will spatially evolve 
to a marginally stable form. Hence, upstream from the shock the reflected 
ion distribution should approach a quasi-plateau. 
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V. ION BEAM TURBULENCE AND SHOCK STRUCTURE 

In Sec. Ill we showed that the experimental mach number vs shock potential 
relationship was in rough agreement with a steady state shock theory involving 
a Maxwellian distribution for incident and reflected ions and the maximal 
trapping electron equation of state. This theory, however, did not include 
the effects of the unstable ion beam turbulence and the turbulent modification 
of the reflected ion distribution function. Since measurements of the electron 
distribution function were not capable of detecting the presence of trapped 
electrons, the question remains as to whether turbulent ion dynamics could 
explain the Mach number-potential relationship even with isothermal electrons. 

An additional question is to what extent the shock structure calculated with 
trapped electrons is modified by ion acoustic turbulence. 

In order to resolve these questions, in this section we derive a differen- 
tial equation for the steady state shock structure which includes the ion beam 
turbulence. We assume that the unstable ion waves are propagating with a 
finite group velocity in the shock frame. Hence, the calculations apply to 
the second type of turbulent shock of Sec. Ill, but probably do not apply to 
the first type of turbulent shock in which the waves appear to have an almost 
zero group speed in the shock frame. We treat the wave dynamics in the quasi- 
linear approximation using WKB spatial eigenfunctions. Rather than solve for 
the turbulent modification of the reflected ion distribution function, we 
assume that the wave turbulence has flattened the distribution function across 
the resonant velocity region and use flux conservation to determine the velocity 
width of the plateau. This assumption yields an upper limit to the turbulent 


modification of the shock structure. 
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Quasilinear Equations 

We assume that in the co-moving frame the shock has a steady, time- 
averaged spatial structure which depends only on the direction x normal to 
the shock front; x = -» is upstream and x = 0 is the first maximum in the 
steady electrostatic potential. We can then write the Vlasov-Poisson equations 
for the steady (fluctuating) particle distribution functions f(;v,x) [f^(v,r)] 
and potential <j> (x) [<j> (_r , t) ] as 


9 f q .7 3f 
a a 8jj>. a 


x 3x 
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f d 3 v 
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(19) 


where q^ = ±e is the electronic charge for ions (+) and electrons (-) , m^ is 

the specific mass, and £ is the sum over ions and electrons. Equation (18) is 

a 

the linearized Vlasov equation. 


Linear Solution 

For low Mach number shocks, the scale length of the shock front is some- 
what larger than the wavelength of the unstable turbulence, which is comparable 
to the Debye length Xq- Therefore to obtain an approximate equation for the 

shock structure we choose WKB solutions to the linear equation in the form 

r x -i 


- l f <*.v> exp 

CO 


k (x')dx' + ik y - iu)t | 
• 7 


( 20 ) 
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(Jj^Cr, t) = £ <J> w (x) exp i k x (x')dx' + ik y - iwt 

i J x ^ 


where cf>^(x) are slowly varying WKB amplitudes, k^(x) is the complex, 

spatially dependent wave number in the x direction, k^ is the y-direction 

wave number, and 0) is the real eigenfrequency . For simplicity, we only 

consider two-dimensional wave propagation. 

By the WKB assumption, 3/3x(f , d> ) is small compared with k (f , d) ) . 

J r w,a v x v w,a’ T a) 

Therefore, we solve Eq. (18) by iteration by writing f = 

J jo oj,a 

where f^ is order k and f^) orc jer 3/3x. After substitution into (18), 
o),a x u,a 


we find 
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where — h q^(x), W^ = . In Eq. (23), the x dependence of 

f a (Wx,Wy,x) has been neglected since it arises solely from the wave turbulence, 

and hence is second order in <p . After substitution of (21), (22), and (23), 

to 

Eq. (19) can be solved to first order in 3/3x to obtain 
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The lowest order in 3/3x solution e(k,u)) = 0 determines the complex wave 

fk 2 3e 


2 2 1/2 

number k = (k^ + k ) ; the next order in 3/3x specifies that 


l*a.l 

x 


is constant, and thus determines the spatial dependence of the WKB amplitude 
2 

|<{> | . If the imaginary parts of k and e are small, then we only need 
Re = const - 


Equation for the Shock Structure 

We now derive an equation for the shock potential <{>, which includes the 
unstable wave turbulence driven by the reflected ions. We first multiply Eq. 
(16) m a v x > integrate over velocity space, and then sum over species to find 
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where e is a unit vector in the x direction. We now eliminate £ q 
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by Eqs. (17) [(19)], substitute (21) for cj>^(_r,t), and perform 
the time average in the last term. If we then write k^(x) as k x = k^ + ik^, 
expand the last term in (25) in order d/dx taking k^ of order d/dx, and take 
the real part of the result, (25) becomes 
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In (26) we have neglected the electron inertia term and retained only the 

electron pressure P ; for simplicity, we have changed notation to write 
2 r 2 2 

k^ E k^, k = (k^ ) + ky . Equation (26) can be integrated with the boundary 
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condition that d/dx = 0 at x = - 00 
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To determine P , we form the momentum moment of the electron equation 

(16), substitute (22) and (23) for f 1 and perform the time average and the 

l , e 

integrals over velocity to find 
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where n = 
e 


2 - 

d v f g . Since the unstable waves have phase speeds on the order 


, 1/2 


of the ion acoustic speed c = [T /M.] and are driven by ion streaming, 

S 6 1 

Im D g << Re D g ; hence, neglecting Ira D g and taking the real part (28) yields 
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If we assume that the electrons are approximately isothermal P g = n e (x)T £ with 
T g = const, (29) can be integrated with respect to x from - 00 to x 
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2 2 2 

For weak turbulence and low Mach number shocks, e d> I /T << 1. Further- 

1 60 1 e 

more, the x' dependence of the integrand in (30) will be dominated by the 


x 


. — f)j _ 

exponential exp(-2 J k^dx"). Hence, we can take n (x’) ^ n(x) outside the 

—CO 

2 

integrand, and expand the exponential to lowest order in | | to find 
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To determine the ion contribution M. 
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f ,2 2 7 

d v v f . 
x i 


in Eq. (27) we form 

the momentum moment of (16), substitute (22) and (23), time average and take 
the real part of the result to find 
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Substituting the quaslinear relation k. = -Im D./(3Ree/3k ) and 

11 X 

integrating (17) over x we find 
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where ru has been written as a function of the potential cf>. 

We can now combine (27), (31), and (33) to obtain the desired equation 
for <j). If we normalize the variables as ip = e^/T^, x -»■ x/Ap, k kA Q , and 
ip = ed> / T , we obtain 
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The local dispersion relation e(k,u>) = 0 was used in reducing the turbulence 
terms. Equation (34) is an implicit equation for ip and would be self- 
consistent if the \p dependence of e(k,w) and k^ were known. If we assume 
that k (x) is only weakly dependent on x, we can replace k by an average 

X X 

•k 

value k^ and perform the x integration to obtain 
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( 35 ) 

where V(ip) is the nonlinear quasipotential. V(^) = 0 occurs at the upstream 
point ip = 0 and at the first maximum of the potential <^(0) . The condition for 
a steady shock solution to (35) is that V(i|>) be negative between ip = 0 and 
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ip = ip(0). If V (\p ) for a given Mach number and temperature ratio T^/T^ is 
negative for all \p , then no steady shock solutions exist. 

Isothermal Electrons 

2 

Bardotti and Segre showed that for T g /T^ = 20 and Maxwellian ions, the 

isothermal electron term e^ forced V(\p) to be negative for Mach numbers (M) 

based on the hot ion linear acoustic speed exceeding 1.07. Recall that, 

experimentally, shocks were produced up to Mach numbers of 1.5 and the Mach 

number-potential maximum relationship disagreed with the predictions of 

2 

Bardotti and Segre. We now investigate whether the turbulent modifications 
to V(40 can bring the isothermal electron theory into better agreement with 
experiment . 

First, note that the explicit turbulence term in (35) makes a positive 
definite contribution to -V(\p), as does . Therefore, the turbulence term 
acts to lower the critical Mach number above which steady shock formation 
is impossible. The possibility remains, however that turbulent modifications 
to the ion distribution function might permit shocks at higher Mach numbers. 

We therefore construct a model for the reflected ions which maximize the 
turbulent modification of the ion distribution function. 

To determine n^(ip) would require solving Eq. (16) for the ion distribution 
f^. However, the convecting wave turbulence resonantly interacts primarily 
with the reflected ions (v < 0) while the incident ions (v^ > 0) are only 
slightly affected by the waves. Therefore, a reasonable approximation for 
the incident ion distribution function f^v^.x) is to assume that f > is just 
the unperturbed upstream distribution. If we assume that f at x = -°° is a 
drifting Maxwellian, then after integration over v^ we can write 
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i 1/2 2' 

f > (V x’ x) = A exp — 2 7 ” U (T* (36) 

a . 

i 

where is the upstream flow speed (or shock speed in the laboratory frame), 
and A is a normalization constant to be determined. 

The instability analysis of the truncated Maxwellian distribution 
indicates that the spatial growth rate is very large, especially in the shock 
front. Therefore, we would expect a very rapid turbulent relaxation of the 
reflected ion distribution f^Cw^jx) to a marginally stable state, which would 
be well-approximated by a plateau in the velocity range where v^ is comparable 
to the unstable wave phase speeds w/k^. In a more exact treatment, resonance 
broadening would enhance the formation of a plateau, and permit diffusion to 
velocities exceeding the specular edge velocity -v^(x) of the unperturbed 
reflected Maxwellian. 

Since ions with v^ = 0 should undergo little diffusion, the reflected and 

incident distributions will be equal at v =0. If we assume that wave 

x 

turbulence flattens f < from v^ = 0 to a new edge velocity -v^(x) and that 
— 

f % 0 for v < -v, (x) , we can write 
< x 1 ’ 

f (W ,x) = f (v = 0,x) = f [W = ecj>(x)] . (37) 

^ X / X A 

To determine v^(x), we can use the conservation of particle flux which is an 
exact consequence of Eq. (16). Steady state requires that the reflected ion 
flux J^(x) at each x be equal to the incident ion flux in the velocity range 
0 £ v^ _< Vq(x), i.e., incident ions which will eventually reflect. J x (x) is 
given by 

, v o <x) 4(0) 

J (x) - v f (W )dv - i- I fju )dW . (38) 

x x > x x M. J_ > x x 

0 1 e<{>(x) 
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The reflected ion flux is 
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The incident (n^) and reflected (n < ) ion densities are 
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If we set the upstream [(j>(x = -«>) = 0] ion density n (0) + n^O) equal to the 
upstream electron density n^(x = -°°) , and substitute (36) into (41), the 
normalization constant A is given by 
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Equations (36) and (41) were substituted into (35) and the critical Mach 


( 42 ) 


number for T /T. = 20 was determined as M = 1.1, or roughly identical to the 
e l 
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Bardotti and Segre result. In this calculation the explicit turbulence term 
in (35) was neglected since the plateau distribution for the reflective ions 
already includes the maximal turbulent modification of the ion distribution 
function. Since this critical Mach number is well below the experimentally 
observed Mach numbers, we conclude that the experimental shocks are inconsistent 
with the theoretical assumption of isothermal electrons even when turbulent 
modification of the ion distribution functions is included. 


Maximally Trapped Electrons 

We now investigate the turbulent modifications to the theoretical shock 
structure predicted by the maximally trapped electron equation of state, Eq. 
(2), We replace the isothermal term e 1 ^ in (35) by P(\p)/P^ from (2) and solve 
for the Mach number as a function of the maximum potential ^ (0) . As a first 
approximation, we assume that the turbulence only slightly modifies the 
reflected ion distribution so that the M vs i^(0) relationship should only be 
slightly shifted from that computed assuming Maxwellian ions. If we expand 

about the value of determined by setting the first two terms in (35) 
equal to zero, the shift in the maximum potential Alp = ip(0) - due to the 
explicit turbulence term in (35) is given by 
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where jjh I = I di I [exp (-2 f k.dx') - 1] is the maximum wave amplitude at 
1 1 max 1 r u ) 1 i l 

— OO 

x = 0 and n > (ip m ) is given by (41). 

For T g /T^ = 20 and Mach numbers 1.1 £ M £ 1.4, the denominator in (43) is 

approximately -0.65. The maximum wave amplitudes observed in the experiments 

are on the order of lib I ol 10 \ From Eq. (24) it can be shown that 

1 1 max 

3Ree/3k is approximately 3Re£/3k n! 2/k [1 + 2(k ^/k^)] nJ 4/k . For the 

X XXX X 

f\j r\j 

oblique unstable ion waves k ^ 1.5, so that we obtain Alp 'v 0.07. Figure 15 
shows the M vs ip(0) relation for both the Maxwellian ions without turbulence 
and the same curve shifted by Aijj = 0.07. Since the A^-shifted curve is in 
poorer agreement with experiments, we conclude that turbulent modification of 
the reflected ion distribution makes a significant change in the shock structure. 

We also calculated the M vs ijj(0) curve for the plateau reflected ion 
distribution derived above, maximally trapped electrons, T g /T^ = 20, and 
neglecting the explicit turbulence term in (35). This result is also shown in 
Figure 15. From the agreement between the theoretical and experimental curves 
we conclude that the electron distribution in the experimental shocks should 
closely approximate the maximally trapped electron distribution and that the 
dominant effect of the wave turbulence on the shock structure is well repre- 
sented by a flattened reflected ion distribution. 

Beam Decay Distance 

Finally, we estimate the spaital extent of the region in which the reflected 
ions interact strongly with the ion acoustic turbulence, which gives a measure 
of the turbulent shock thickness. Separating the reflected ion contribution 
from Eq. (32) and neglecting all contributions proportional to the reflected 
ion density n < except the kinetic term, we have 



-31- 


d_ 

dx 


, 2, 3, ”<"A > 

(n <M+ v x ) -v 

J5 



CO 


k k 2 k. 
x 1 

4tt 


T 

3Ree e 1 ,2 

3k 2 * co max 

x e 


(44) 


where is an effective beam decay distance and <p^ has been replaced by its 
maximum value. In dimensionless units (44) becomes 
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For typical experimental parameters n < /n e ^ 0.05, <v^ > 'v 4c g , k^ 'V 0.1, 

I^J 2 £ 10' 2 , and k = 1.5, we find L^/Ap ^20, in rough agreement with experi- 
mentally observed values of the spatial extent of the wave turbulence region. 
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VI . SUMMARY 

Electrostatic ion acoustic shock experiments in the University of 
California, Los Angeles, double plasma device have revealed three distinct types 
of shock structures. 1) A laminar shock with a trailing wave train of potential 
oscillations in which few ions are reflected from the shock potential jump; 

2) a weakly turbulent shock in which the reflected ion-driven ion acoustic 
turbulence is confined ahead of the major jump in shock potential; 3) a 
more strongly turbulent shock in which the ion acoustic turbulence propagates 
throughout the shock front and has amplitudes comparable to the steady shock 
potential. In the first two types of shocks, the potential jump is sharp 
occurring on scale lengths of a few Debye lengths; for the third type of shock, 
the shock transition is considerably broader, and occurs on scale lengths on 

the order of the reflected ion decay distance, some tens of Debye lengths. 

The unstable ion acoustic waves are believed to propagate at large angles to 
the shock normal. Figure 16 presents a qualitative summary of the shock in 
terms of the parameter range of the experiments. The boundaries between the 
small and large amplitude turbulent shocks are not precise. Also recall that 

the shock evolves from one state to another as it propagates away from the 

formation region. 

For T^/T^ in the range 10 to 40, shocks with Mach numbers between 1.05 

and 1.5 have been observed. These Mach numbers generally exceed the maximum 

2 

Mach number calculated by Bardotti and Segre for isothermal electrons and 
Maxwellian ions. Steady state shock calculations assuming trapped electrons^ 
and turbulently flattened reflected ions predict a Mach number vs maximum shock 
potential relationship which agrees well with experiment. This agreement 
suggests that in the double plasma shock experiments the electrons are indeed trapped 
in the steady shock potential. Electrons could become trapped between the 
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negatively biased excitation grid and the shock front, and execute several 
bounces in this effective potential well before colliding with the side chamber 
walls. Hence, the existence of a flat-topped distribution of trapped electrons 
in these experiments is probably reasonable. Indeed such a trapped-electron 
distribution was recently observed (in another UCLA UP device) in the potential 

f^l % 

maxima of current-driven ion acoustic waves of large amplitude I— ^ 102 

' o 

The difference between the two types of turbulent shocks is not currently 

understood. The weakly turbulent shocks appear to require that the unstable 

ion acoustic waves be close to group standing in the shock front. The linear 

instability theory of a Maxwellian ion beam plasma indicates that the group 

standing condition might be satisfied, although a parametric investigation of 

the group standing condition is prohibited by our imprecise knowledge of the 

reflected ion distribution functions and the angular distribution of the 

unstable waves. If the waves have an approximately zero group velocity in the 

shock frame, then wave saturation becomes an interesting theoretical and 

experimental problem. Since in steady state these waves would be continuously 

driven unstable by the reflected ions, wave saturation could occur by mode 

coupling to non-group standing waves, as in the early perpendicular magnetosonic 

23 

whistler shock theory of Camac et al. Alternatively, if the waves propagated 
slightly faster than the shock front, the waves might simply damp out ahead of 
the shock. A more detailed experimental and theoretical investigation is 
needed to clarify these points. 

In the strongly turbulent shocks, ions are no longer reflected by a smooth 
laminar potential barrier but by a spatially and temporally turbulent shock 
potential. Hence, the essentially laminar theory, which was used in this paper, 
is probably inadequate to accurately describe the shock structure, and effects 
such as ion trapping in the shock front must be considered. At what point the 
essentially fluid shock description with superimposed wave turbulence goes over 
into the fully turbulent kinetic shock description remains to be resolved. 
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Fi£ 1 


Fig. 2 


Fig. 3 


Fig. 4: 


Fig. 5: 


Figure Captions 

Laminar shock profile of electron density AN/N vs time as a function 

of distance x from the separation grid of the DP device. The foot 

structure and reflected ions can be seen ahead of the sharp jump in 

electron density. Downstream of the shock density jump there is a 

trailing wave train of density oscillations. M = 1.13, Tg/T^ = 13, 

V L /C = 1.6. 
b s 

Mach number M vs maximum shock potential etj>/Tg. The crosses are 
experimental points for shocks with 5 Tg/T^ — 50. The solid curve 
is the average of the experimental points M = 1.0 + 0.6 e<j>/Tg. The 
lower dashed curve is computed from the isothermal-cold-ion Sagdeev- 
Moiseev theory. The upper dashed curve is computed using the trapped 
electron equation of state with cold ions. 

Mach number M vs maximum shock potential ecJ>/Tg calculated using the 

trapped electron equation of state and a Maxwellian ion distribution 

for the incident ions and a truncated Maxwellian for the reflected ions. 

The density of ions reflected from the shock front n, /n~ vs the 

b u 

maximum shock potential ed>/ T for T /T. = 10 and 30. The reflected 

e e x 

ion density predicted from a truncated Maxwellian reflected ion distri- 
bution and the experimental M vs the maximum etjj/T^ is shown as the 
solid line. The reflected ion density predicted by the Sagdeev-Moiseev 
M vs the maximum e^/T^ is shown as the dashed line. 

The reflected ion density mean velocity V^/Cg vs the maximum shock 
potential ec|>/Tg for 5 £ Tg/T^ _< 50. The solid curves were calculated 
assuming a truncated Maxwellian for the reflected ion distribution and 
the experimental average Mach number maximum shock potential relation. 



-38- 


Fig. 6: The electron density AN/N vs time as a function of distance x from 

the separation grid for a small amplitude turbulent shock. The 

turbulence appears as the blurred region. M = 1.2, V^/C = 1.8, 

T /T. = 10, n,/n_ = 0.08. 
e l bU 

Fig. 7: The electron density AN/N vs time as a function of distance x from 

the separation grid for a large amplitude turbulent shock. M = 1.2, 

V v /C = 2.0, T /T. - 16, n, /n_ = 0.05. 
b s e l b 0 

Fig. 8: Individual large amplitude turbulent shock electron density profiles 

vs time. Each oscilloscope trace consists of three shock profiles 
taken at the same distance from the DP separation grid. The shock 
parameters are the same as Fig. 7. 

Fig. 9: Continuation of Fig. 8 for further distances from the separation grid. 

Fig. 10: The density ratio n(x) defined as the local shock electron density 
jump normalized to the initial shock electron density jump at x = 2 
vs distance x. 


Fig. 11: The angle 0 to the direction of the ion beam at which unstable ion 

acoustic waves are detected vs the ion beam energy etJ^/T • The 
cross-hatched region represents a 3 db width in the power spectrum. 

The two solid curves represent approximate theoretical stability limits. 


The maximum spatial growth increment k^/k^ for ion acoustic waves as 
a function of beam density n^/n^ for various ij T^ ratios. The beam 
velocity is 1.0 C g . 

The group velocity V /C in the direction of the ion beam vs frequency 

g s 

ui/uipi* The lower curve is calculated neglecting the beam contribution 
to the group velocity; the upper curve includes the beam contribution. 


Fig. 14: The spatial growth increment Im k^ vs potential e<j>(x)/T e in the shock 


front. 
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Fig. 15: Mach number M vs the maximum e^/T^. The dashed curve is calculated 

with the electron equation of state and a truncated Maxwellian for 

the reflected ions; T /T. = 20. The lowest curve is the dashed 

e 1 

curve shifted by A(e<j>/T e ) = Alp = 0.07. The average experimental M 
vs efJi/T^ relationship is shown as the top curve. The intermediate 
curve is calculated assuming trapped electrons and turbulently 
flattened reflected ions. 

Fig. 16: Schematic summary of experimental shock results presented in a 

maximum shock potential vs T^/T^ parameter space. The boundaries 


between the various types of observed shocks are not precise. 
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